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1  1.  INTRODUCTION 

2  The  purpose  of  this  paper  is  to  build  on  the  hierarchical  framework  and  concepts  described  in 

3  a  companion  paper  to  develop  a  methodology  for  the  Indicator  of  Risk  of  Water  Contamination 

4  (IROWC).1    The  methodology  is  proposed  as  a  preliminary  model  only.    It  requires  further 

5  development  and  continued  testing  to  authenticate  its  utility  among  hierarchical  levels,  and  at 

6  various  temporal  and  spatial  scales. 

7  This  discussion  document  describes  the  initial  development  and  demonstration  of  IROWC,  and 

8  suggests  areas  of  further  research,  testing  and  implementation.    The  paper  first  briefly  reviews 

9  a  number  of  risk-based  approaches  to  the  analysis  and  assessment  of  water  contamination  by 

10  agriculture,  based  on  previous  research.    It  then  proposes  a  methodology  for  IROWC  using  the 

1 1  budget  approach.    The  proposed  methodology  is  demonstrated  utilizing  preliminary  data  and. 

12  finally,  suggestions  for  detailed  implementation  are  presented. 

13  2.  RISK-BASED  APPROACHES 

1 4  The  potential  of  nitrates  or  pesticides  to  contaminate  water  is  frequently  determined  by  screening 

15  or  simulation  models.     Screening  models  assess  and  rank  the  leaching  potential  of  inputs. 

16  generally  on  the  basis  of  chemical  properties.    Simulation  models  incorporate  various  physical. 

1 7  chemical  and  biological  attributes  of  the  input  and  site  to  predict  the  amount  of  a  contaminant 

18  which  is  available  for  contamination.    Both  models  are  apparent  in  risk-based  approaches. 

19  This   section   identifies  and  describes  previous  risk-based  approaches  to  the  analysis  and 

20  assessment  of  nitrate  and  pesticide  contamination  of  surface  and  ground  water.   It  classifies  the 

21  approaches  into  general  categories  and  assesses  their  utility  for  IROWC  using  the  desirable 

22  characteristics  of  IROWC  identified  in  the  companion  paper. 

23  2.1  RISK  ANALYSIS  USING  SPATIAL  INDICATORS 

24  Spatial  indicators  generally  map  spatial  attributes  of  IROWC  components  (e.g.,  precipitation,  soil 

25  leaching  potential,  crop  type)  to  identify  areas  at  risk  of  contamination.    GIS  is  a  common  tool 

26  for  analysis  and  display. 


MacDonald,  K.B.  and  H.  Spaling.  1995.  Indicator  of  Risk  of  Water  Contamination:  Concepts 
and  Principles.    Working  Paper.  Draft  of  21  March  1995,  Agriculture  and  Agri-Food  Canada,  Ontario  Land 
29  Resource  Unit,  Centre  for  Land  and  Biological  Resources  Research,  Guelph,  Ontario. 
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1  2.1.1  Vulnerability  Mapping 

2  Vulnerability  maps  identify  and  rate  areas  according  to  contamination  potential.    Vulnerability 

3  is  based  on  the  spatial  association  of  attributes  which  affect  contamination.    Spatial  analysis  is 

4  usually  integrated  with  a  contaminant  screening  or  modelling  procedure. 

5  A  USDA  study  (Kellogg  et  al.  1992)  developed  vulnerability  indexes  for  pesticide  and  nitrate 

6  contamination  of  shallow  ground  water.    The  Ground  Water  Vulnerability  Index  for  Pesticides 
(GWVIP)  is  a  function  of  soil-pesticide  leaching,  precipitation,  and  chemical  use.    The  nitrate 

8  index  (GWVIN)  is  based  on  precipitation  and  excess  nitrogen  fertilizer  (i.e.,  applied  minus  crop 

9  uptake  and  removal).    Index  scores  are  grouped  into  relative  categories  of  risk  (high,  moderate, 

10  low).    Data  are  displayed  at  the  national  scale  with  county  level  resolution. 

1 1  McRae  (1991)  also  uses  a  vulnerability  mapping  procedure  as  part  of  a  study  to  identify  areas 

12  at  risk  of  ground  water  contamination  in  Canada.    It  is  less  comprehensive  than  the  USDA 

13  approach.      The   focus   is   on  pesticide   leaching   potential,   using   only    soil   and   landform 

1 4  characteristics  derived  from  soil  polygons  on  the  Generalized  Soil  Landscape  Maps  (1:1 ,000,000). 

15  Areas  vulnerable  to  pesticide  leaching  are  characterized  by  sandy  or  sandy  loam  classes  of  soil 

1 6  texture,  a  slope  gradient  of  0-9%,  and  undulating,  hummocky,  level,  or  knoll  and  kettle  classes 

1 7  of  surface  formation.  McRae  ( 1 99 1 )  suggests  improving  the  procedure  by  including  ground  water 

18  recharge  data. 

19  The  above  studies  generalize  on  the  basis  of  soil  type  (e.g.,  sandy  soils  are  equated  with  high 

20  leaching  potential),  but  caution  is  needed  since  leaching  may  be  higher  in  some  instances  under 

21  other  soils  and  certain  management  conditions  (e.g.,  clay  soils  under  ridge  or  no-tillage  where 

22  preferential  flow  is  greater). 

23  In  another  application,  Khan  and  Liang  (1989)  use  a  GIS-based  approach  to  construct  maps 

24  showing  the  likelihood  of  pesticide  contamination  of  ground  water  for  Hawaii  (1500  km2).    A 

25  "contamination  likelihood  map"  for  various  pesticides  is  based  on  an  attenuation  factor  (AF) 

26  derived  from  the  percent  of  applied  chemical  which  leaches  below  the  surface  soil  layers.  A  GIS 

27  is  used  to  delineate  the  land  area  falling  into  different  AF  rating  groups,  and  to  calculate 

28  differences  in  contamination  potential  of  the  same  pesticide  for  different  crops. 

29  These  studies  show  that  vulnerability  mapping  is  capable  of  integrating  IROWC  components  by 

30  determining  their  spatial  association.  The  geographic  scale  and  resolution  of  these  studies  suggest 

31  that  the  procedure  is  most  useful  at  levels  4-3  or  higher.    The  analysis  can  incorporate  a  wide 

32  range  of  spatial  attributes,  although  data  requirements  are  extensive.   The  product  (map)  is  easily 

33  communicated  and  interpreted.    Applications  focus  on  the  risk  of  ground  water  contamination, 

34  but  the  procedure  could  be  adapted  for  surface  water. 
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1  2.1.2  Risk  Based  on  Buffer  Strip  Characterization 

2  This  approach  analyzes  the  width  of  an  area  adjacent  to  surface  waters  (riparian  zone)  which 

3  acts  as  a  protective  strip  to  delay,  assimilate  or  filter  contaminants  in  runoff  before  they  enter 

4  a  lotic  system. 

5  A  buffer  strip  is  a  permanently  vegetated  area  of  land  situated  between  the  contaminant  source 

6  areas  (i.e.,  agricultural  fields)  and  the  receiving  waters.     Buffer  strips  have  the  potential  to 

7  reduce  contamination  from  diffuse  sources  through  the  retention  of  nitrogen,  phosphorous  and 

8  pesticides  in  runoff.     Several  reviews  have  assessed  the  effectiveness  of  buffer  strips  as  a 

9  mechanism   to  reduce  surface  water  contamination   from  agriculture   (Muscutt   et  al.    1993. 

10  Osborne  and  Kovacic  1993.  Phillips  1989). 

11  Effectiveness  (i.e.,  contaminant  retention  potential)  is  directly  related  to  width  of  the  buffer 

12  strip,  and  the  strip's  environmental  attributes  and  management  characteristics.     Width  can  be 

13  determined  by  i)  defining  a  constant  distance  (i.e.,  fixed  setback  from  all  surface  waters),  ii) 

14  establishing  a  minimum-variable  distance  (i.e.,  beyond  a  minimum,  width  varies  depending  on 

15  environmental    conditions    and    land    use),    or    iii)    using    a    variable    distance    based    on 

16  environmental  and  land  use  factors.    Variable  buffer  widths  have  been  calculated  for  streams 

17  in  a  rural  county  of  N.  Carolina  (Xiang  1993).    This  study  incorporates  a  pollutant  detention 

18  time  model  into  a  GIS  framework  to  analyze  effectiveness  of  buffer  strip  widths  based  on  soil 

19  hydrology,  land  cover  and  topography. 

20  Buffer  strips  are  most  suitable  as  indicators  of  risk  for  surface  water.     Risk  analysis  •  and 

2 1  assessment  could  be  based  on  deviation  from  an  optimum  buffer  width.   Various  contaminants 

22  could  be  integrated  into  one  indicator  by  identifying  the  contaminant  with  the  widest  buffer 
requirement.   Like  vulnerability  mapping,  IROWC  based  on  buffer  strips  is  conducive  to  spatial 

24  presentation,  which  is  easily  communicated  and  interpreted.    A  need  for  comprehensive  data 

25  bases  suggests  that  application  is  most  appropriate  at  level  3,  and  perhaps  4.    A  significant 

26  shortcoming  is  that  subsurface  drainage  (tile  drains)  is  not  accounted  for  and  could  lead  to  an 

27  error  in  the  calculation  of  risk  for  high  rainfall  areas  (e.g.,  eastern  Canada). 

28  2.2  BUDGET  APPROACHES  TO  RISK 

Budget  approaches  use  input-output  analyses  to  determine  the  amount  of  a  plant  nutrient  or  a 

30  pesticide  which  is  in  excess  of  crop  requirements  and  may  be  available  to  contaminate  water. 

31  Barry  et  al.  (1993)  use  this  approach  to  predict  nitrogen  balances  for  various  crop  rotations  in 
southwestern  and  western  Ontario,  and  two  farming  systems  (cash  crop  and  dairy)  in  western 

33  Ontario.   The  quantity  of  surplus  nitrogen  is  adjusted  by  an  average  ground  water  recharge  rate 

to  estimate  the  potential  NO--N  concentration  in  shallow  ground  water  (assumes  an  average 
ground  water  recharge  of  160  mm  yr'1  for  southern  Ontario  based  on  a  proportion  (51.8%)  of 
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1  the  long  term  mean  annual  discharge  of  ten  agricultural  watersheds).    For  example,  estimated 

2  nitrate    concentrations    in   ground    water    under    a   corn-soybean-winter    wheat    rotation    in 

3  southwestern  Ontario  was  8.5  mg  N03-N  L"1  and  22.4  mg  N03-N  L"1  in  western  Ontario  for 

4  the  same  cropping  rotation,  and  58.4  mg  L"1  for  a  dairy  farm. 

5  A  major  initiative  to  determine  annual  nutrient  budgets  for  N,  P  and  K  is  underway  for  several 

6  locations  in  the  lower  Fraser  Valley  of  British  Columbia  (B.  Zebarth,  personal  communication, 

7  1995).    An  annual  nutrient  balance  (applied  surplus  or  deficit)  is  based  on  type  and  number 

8  livestock,  manure  management,  inorganic  fertilizer  applications,  crop  nutrient  requirements,  and 

9  atmospheric  input.   Estimates  of  soil  release  pathways  include  surface  and  ground  water.    Data 

10  on  inputs,  pathways  and  outputs  are  differentiated  on  a  monthly  basis. 

1 1  Budget  approaches   are  also  represented   in  calculations   of  mass  balance  to  determine  the 

12  dissipation   of  pesticides.      For  example,  Frank  et  al.  (1991a,    1991b)  have  calculated   the 

13  proportion   of  applied   atrazine,   cyanazine   and  metolachlor   which   is  dissipated   by  runoff. 

14  subsurface  tile  drainage,  leaching  to  shallow  ground  water,  and  degradation  in  clay  loam  soils 

15  of  Ontario  (field  scale). 

1 6  Risk  is  related  to  the  quantity  of  a  potential  contaminant  present  in  excess  of  a  desirable  or 

17  tolerable  quantity.     Risk  is  increased  if  a  budget  balance  shows  an  excess  amount.     Unlike 

18  indicators  which  focus  on  inputs  or  outputs  (e.g.,  application  rates  per  hectare,  yield),  a  budget 

19  approach  to  IROWC  considers  the  excess  amount  available  for  contamination2.    The  above 

20  studies  by  Barry  et  al.  (1993)  and  Frank  et  al.  (1991a,  1992b)  apply  a  budget  approach  at  the 

21  farm  and  field  scales  (levels  1  and  2),  but  it  is  potentially  applicable  at  higher  levels,  as 

22  demonstrated    in  the  nutrient  budget   calculations    in  British   Columbia   where   study   areas 

23  represent  an  aggregation  of  farms.    Also,  crude  nutrient  budgets  have  been  estimated  at  the 

24  national  level  (see  Table  9  in  the  concept  paper).    Factors  which  alter  the  budget  balance  are 

25  readily  incorporated  (e.g.,  increasing  use  of  conservation  tillage).    Notwithstanding  the  earlier 

26  reference  to  spatial  indicators,  a  budget  approach  can  be  integrated  with  spatial  analysis  to 

27  differentiate  areas  at  risk  of  surface  or  ground  water  contamination  from  nitrogen  or  pesticides. 

28  2.3  RISK  AS  PROBABILITY  OF  EXCEEDENCE 

29  This  approach  calculates  an  estimate  of  pesticide  contamination  in  water,  and  expresses  risk 

30  as  a  probability  of  exceeding  an  established  or  predefined  water  quality  standard  over  a  defined 

3 1  period  of  time. 


32  2  In  areas  with  an  annual  deficit  balance,  unusual  events  may  still  contribute  to  water 

contamination    (e.g.,  high   rainfall   event   (>50mm)   on  tile  drained    soils),   but  the   relative 
contribution  is  likely  to  be  considerably  less  than  that  of  a  surplus  balance.  j 
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1  Varshney    et  al.   (1993)   used   this  approach   in  a  risk-based    evaluation    of  ground   water 

2  contamination  by  atrazine,  simazine  and  alachlor  (southeast  Iowa).     The  simulation  model 

3  RUSTIC  provides  estimates  of  concentrations  of  each  herbicide  within  the  aquifer  (27-yearsi. 

4  from  which  cumulative  frequency  distributions  of  herbicide  concentrations  at  various  locations 

5  are  calculated.    The  potential  risk  (R)  associated  with  each  herbicide  is  defined  as: 

6  R  =  Pr(OCs) 

where  Pr  =  probability  of  exceedence,  C  =  predicted  daily  pesticide  concentration  in  ground 

8  water  (ppb),  and  Cs  =  maximum  contaminant  level  (ppb)  (EPA  standards  for  each  herbicide). 

9  For  example,  simazine  exceeded  Cs  35%  of  the  time. 

10  This  approach  is  attractive   for  IROWC   because   it  explicitly   incorporates  a  water  quality 

1 1  standard,  although  there  is  no  guidance  for  its  selection.     The  use  of  a  simulation  model 

12  (RUSTIC)  to  estimate  pesticide  contamination  is  dependent  on  a  long  historical  record.    For 

13  situations  of  insufficient  data,  alternative  probability  distributions  for  estimating  the  probability 

14  of  exceedence  of  contaminant  concentrations  have  been  suggested  by  McBean  and  Rovers 

15  (1992).   Even  if  temporal  limitations  are  overcome,  the  geographic  application  of  this  approach 

16  is  likely  limited  to  level  2,  and  maybe  3.   Expressing  risk  as  a  probability  of  exceeding  a  water 

17  quality  standard  is  likely  to  be  understood  by  non-technical  «people. 

18  2.4  RISK  IN  RATING  INDICES/INDEXES 

19  Rating  approaches  rank  contaminants  according  to  their  leaching  potential.    Ranking  schemes 

20  for  nitrates  are  based  mainly  on  soil  properties,  whereas  those  for  pesticides  consider  the 

21  physical    and   chemical    properties   of  the   material,    as   well   as   other   soil,    landscape    and 

22  environmental  characteristics.    Land  use  and  management  are  frequently  held  constant. 

23  Khakural   and  Robert  (1993)   used   LEACHM-N   as  a  screening   tool  to  develop   soil  NO, 

24  leaching  potential  (NLP)  ratings  from  soil  survey  information  and  representative  weather  data. 

25  Seasonally  accumulated  values  of  leached  N03  (>1.5  m)  for  a  normal  precipitation  year  were 

26  used  as  NLP  indices.    The  ratings  were  added  to  the  Soil  Survey  Information  System  (SSIS) 

27  to  create  N03  leaching  potential  rating  maps  for  three  Minnesota  counties. 

28  McRae  (1991)  developed  a  pesticide  rating  scheme  for  Canada  which  ranked  pesticides  based 

29  on  three  factors:  leaching  potential,  volume  of  use,  and  overall  regulatory  concern.    Pesticides 

30  are  ranked  according  to  their  leaching  potential  using  Gustafson's  technique  (based  on  partition 

31  coefficient  between  organic  carbon  and  water,  and  half-life  of  the  pesticide  in  soil).    Rankings 

32  by  volume  of  use  refer  to  kg/annum   of  active  ingredient.      Ranking   according   to  overall 

33  regulatory  concern  reflects  a  pesticide's  status  in  the  pesticide  re-evaluation  program.     The 

34  scheme  does  not  integrate  the  three  rankings,  but  prompts  the  user  to  consider  each  ranking  in 

35  a  sequential  manner. 
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1  Other  pesticide  rating  schemes  also  have  been  developed.     For  example,   Warner  (1985) 

2  proposed  an  environmental  risk  index  which  assigns  a  relative  value  to  pesticides  based  on  a 

3  tradeoff  between  its  likelihood  to  contaminate  (depending  on  persistence  and  mobility),  and  its 

4  toxicity.    The  index  differentiates  risk  for  surface  and  ground  water,  and  is  applied  to  various 

5  agricultural  crops  in  Suffolk  County,  NY. 

6  Rating  approaches  primarily  focus  on  properties  of  contaminants  and/or  soils.   For  this  reason, 

7  their  utility  for  IROWC  is  limited,  unless  integrated  into,  or  accompanied  by,  approaches  (e.g., 

8  spatial    analysis)    which    consider    other    IROWC    components    (e.g.,    climate,    land    use, 

9  management). 

10  2.5  RISK  AND  EXPERT  SYSTEMS 

1 1  Expert  systems  are  computer  models  which  help  non-experts  to  find  solutions  to  complex 

1 2  problems.   The  models  require  not  only  quantitative  data  but  also  qualitative  information  in  the 

13  form  of  decision  rules,  interpretation  and  judgement. 

14  Crowe  and  Mutch  (1994)  have  developed  an  expert  systems  approach  (EXPRES)  to  evaluate 

15  the  potential  of  a  pesticide  to  contaminant  ground  water  (see  also  Crowe  1994).    Depending 

16  on  study  objectives,  available  data  and  time  constraints,  EXPRES  helps  the  user  to  select  a 

17  pesticide  simulation  model,  prepare  an  input  data  set,  carry  out  the  assessment,  and  interpret 

18  the  results.     The  model  provides  a  relative  assessment  of  the  potential  of  each  pesticide  to 

19  contaminate  ground  water,  including  a  quantitative  prediction  of  the  pesticide  distribution  and 

20  migration  rates  (time  and  depth).    The  model  is  applied  to  a  "typical"  field  of  potato,  irrigated 

21  sugar  beet,  and  corn  in  PEI.  Alberta  and  Ontario,  respectively.      (The  model  contains  22 

22  agricultural  regions  in  Canada). 

23  Although   EXPRES    is  designed   primarily   as  a  screening   tool   for  decision   making   (e.g., 

24  pesticide  regulation),  it  has  significant  potential  as  an  approach  to  IROWC.     Presumably  a 

25  user's  objectives  can  be  defined  in  terms  of  risk  assessment.    The  main  components  of  IROWC 

26  are  incorporated   in  the  model.     However,  the  model's  field-specific  focus  is  likely  to  limit 

27  application  to  levels  1  and  2.  although  aggregation  to  level  3  and  possibly  4  may  be  possible 

28  for  homogeneous  conditions  at  these  higher  levels.     The  model  requires  detailed  data  sets. 

29  While  representative  data  sets  are  available  with  the  model,  ideally  these  should  be  generated 

30  for  each  application. 

31  2.6  OVERALL  ASSESSMENT 

32  The  approaches  reviewed  here  are  relevant  to  the  needs  of  IROWC  and  generally  define  some 

33  potential  methodological  approaches.    However,  with  the  exception  of  the  budgeting  approach, 

34  they  fail  to  incorporate  some  characteristics  of  IROWC  identified  as  desirable  in  the  concept 

35  paper.    Specifically,  most  approaches  provide  a  static  assessment  which  does  not  incorporate 
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1  specific  farm  management  and  production  practices  and  does  not  readily  allow  for  temporal 

2  analysis. 

3  3.  A  PROPOSED  METHODOLOGY  FOR  IROWC 

4  Risk  of  water  contamination  is  based  on  a  concentration  of  a  potential  contaminant  in  water 

5  over  time.    Quite  apart  from  the  criteria  used  to  define  contamination,  concentration  depends 

6  on  two  variables: 

7  1 .     the  quantity  of  a  potential  contaminant,  and 

8  2.     the  amount  of  water  available  to  dilute  it. 

9  The  amount  of  water  available  to  dilute  a  potential  contaminant   is  determined  by  water 

10  partitioning.    At  any  site,  total  precipitation  is  partitioned  into  different  forms.    For  example. 

1 1  during  the  growing  season,  a  certain  proportion  of  water  is  consumed  by  crops,  and  another 

12  is  "lost"  by  evaporation  from  the  soil  surface.     The  remainder,  or  excess  water  above  soil 

13  storage  capacity,  moves  off  the  site  to  recharge  ground  and  surface  water  supplies. 

14  3.1  CONSIDERATIONS   FOR  DEVELOPMENT  OF  THE  METHODOLOGY 

15  3.1.1  Hydrologic  Considerations 

16  In  the  proposed  methodology,  excess  water  plays  a  fundamental  role  in  determining  the  level 

1 7  of  risk  of  contamination.    A  water  balance  or  budget  approach  can  be  used  to  determine  the 

18  quantity  of  excess  water  which  is  available  to  transport  any  residual  amount  of  a  contaminant. 

19  The  basis  for  this  approach  is  the  terrestrial  hydrological  cycle.    In  this  cycle,  the  main  input 

20  (precipitation)  is  transformed  through  various  processes  (infiltration,  plant  uptake,  runoff)  into 

21  output  (evapotranspiration). 

22  Ideally,  the  water  balance  of  an  area  is  represented  by  the  difference  between  precipitation  (P) 

23  and  actual  evapotranspiration  (AE).  However,  data  for  AE  is  frequently  lacking.   It  is  possible 

24  to  determine  a  crude  water  balance  at  any  location  based  on  annual  or  seasonal  values  of  P  and 

25  potential  evapotranspiration  (PE).    When  P  exceeds  PE,  surplus  moisture  is  assumed  to  exist 

26  and  water  accumulates   in  the  soil  until  field  capacity   is  reached.     Any  further  excess   is 

27  available  as  surface  water  (i.e.,  runoff,  subsurface  tile  flow),  or  ground  water  recharge.    When 

28  P<PE,  a  water  deficit  occurs,  and  there  is  no  surplus  moisture  available. 

29  A  water  balance  approach  can  contribute  to  the  proposed  methodology  in  three  ways.    First. 

30  data  on  excess  water  can  be  distributed  over  time  throughout  the  year.    Seasonal  or  monthly 

31  values   of  water   surplus   or  deficit  can   be  calculated,    and  analyzed   against   the  temporal 

32  distribution  of  excess  nitrogen  or  pesticide  residue.    Second,  the  variation  between  years  can 

33  be  assessed.    Third,  the  water  balance  approach  partitions  output  flow  into  surface  (e.g.,  runoff. 

34  seepage)  and  ground  water  (e.g.,  infiltration,  deep  percolation)   components.      Data  on  the 
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1  proportion    of  flow  through   each   pathway    is   needed   to  analyze   and  assess   the   risk  of 

2  contamination  according  to  destination. 

3  Agriculture  is  directly  dependent  on  the  hydrological  cycle  (e.g.,  soil  moisture  for  crops  and 

4  drinking  water  for  livestock),  but  agriculture  also  alters  the  cycle,  particularly  the  transfer 

5  mechanisms  and  pathways  of  flow.   For  example,  tillage  practices  affect  infiltration  rates,  crop 

6  type  and  planting  pattern  influences  runoff,  and  land  drainage  impacts  percolation  and  lateral 

7  flow.     These  impacts  become  increasingly   important  at  lower  spatial  and  temporal  scales. 

8  Thus,  calculation  of  excess  water  at  these  scales  needs  to  account  for  the  effects  of  land 

9  management  on  the  water  balance  of  an  agricultural  area. 

10  While  the  pathways  of  the  hydrologic  cycle  are  reasonably  well  understood,  good  quantitative 

1 1  estimates  of  the  magnitude  of  flows  in  specific  components  are  very  difficult. 

12  3.1.2  Nutrient  Considerations 

13  Nitrogen  cycles  in  intensive,  arable  agriculture  are  generally  characterized  by  the  addition  of 

14  nitrogen  via  fertilizer  and/or  animal  manures,  and  the  removal  of  nitrogen  in  the  harvested  crop 

15  and  cover  crops.     Other  inputs  include  atmospheric  deposition  and  nitrogen  fixation,  while 

16  losses  occur  through  leaching,  denitrification  and  immobilization.     These  components  of  the 

17  nitrogen  cycle,  and  their  proportionate  amounts,  are  shown  for  wheat  production  in  Figure  1. 

18  Nitrogen  in  the  form  of  fertilizer  and  manures  often  constitutes  a  large  fraction  of  the  total 

19  nitrogen  input.    Generally,  crops  respond  positively  to  increased  fertilizer  and  manure  inputs, 

20  but  this  response  also  increases   nitrogen  removed  in  the  harvested   material.     Losses  via 

21  leaching  are  also  usually  higher  with  increasing   nitrogen  input.     As  fertilizer  nitrogen  is 

22  increased,  the  efficiency  of  nitrogen  use  is  decreased  (i.e.,  percentage  of  total  fertilizer  nitrogen 

23  taken  up  by  the  plant  declines).      Nitrogen   use  efficiencies   for  crops  have  been  broadly 

24  estimated  at  60-50%  (e.g.,  Briggs  and  Courtney  1989,  Juergens-Gschwind   1989).   This  means 

25  that  a  considerable  quantity  of  nitrogen  not  taken  up  by  plants  is  potentially  available  for 

26  leaching  or  other  dissipation  pathways  of  the  nitrogen  cycle. 

27  3.1.3  Pesticide  Considerations 

28  An  ideal  approach  to  analyzing  the  risk  of  pesticide  contamination  first  determines  the  quantity 

29  of  a  pesticide,  and  its  chemical  derivatives,  in  the  soil,  and  then  ascertains  the  susceptibility 

30  of  these    chemicals    to    various    transport    and    transformation    processes    (Mackay    1992). 

31  Following  application,  a  pesticide  is  partitioned  into  various  phases  (e.g.,  vapour  in  soil  air, 

32  sorbed   by  organic   matter,   solution    in   soil   water).      Of  particular   concern    is  the  amount 

33  partitioned  in  mobile  soil  water,  which  represents  the  quantity  potentially  available  for  leaching 

34  or  surface  runoff.  The  proportion  directed  toward  leaching  can  be  relatively  high  for  herbicides 
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1  such  as  atrazine  (Table  1).  The  risk  of  surface  transport  (in  solution  or  attached  to  eroding  soil 

2  particles)  is  especially  high  if  rainfall  closely  follows  pesticide  application. 

3  Quantitative  models  of  chemical  fate  in  soil  are  available  to  estimate  the  amount  of  a  pesticide 

4  partitioned  in  soil  water  (e.g.,  PRZM,  LEACHM).  These  models  are  generally  site-specific  and 

5  applicable  only  at  lower  hierarchical  levels 

6  (levels    3-1).      Their   data   demands    are 

7  extensive  Table  1            Fate  of  Atrazine  in  Soil 

8  At     higher      levels      (e.g.,      levels      5-3),         Loss  By:  % 

9  information  is  available  on  pesticide  inputs  reaction                     31.4 

10  (e.g.,  recommended  application  rates  by  crop,  leaching                    68.0 

11  surveys  of  pesticide  use),  but  understanding        '     evaporation  to  air      0.6 

12  of  and  data  on  partitioning,  transformation  Source:  data  from  Macka>'  <1991) 

13  and  transport  are  lacking.     For  this  reason,     

14  the   relative    indicator   of  risk   of  pesticide 

15  contamination  at  level  5  considers  pesticide  input  data  only. 

16  3.1.4  Key  Questions 

17  Based  on  the  methodological  considerations  discussed  above,  three  questions  for  investigation 

18  can  be  posed: 

19  1.     is  there  sufficient  excess  water  on  average  to  move  contaminants  (vertically  or 

20  laterally)  to  water  sources  (i.e.,  is  there  a  risk?) 

21  This  question  in  reality  deals  not  only  with  the  amount  of  excess  water,  but  also  the 

22  proximity  of  the  rooting  zone  to  water  resources   (surface   or  ground).      It  can  be 

23  rephrased  as  "what  is  the  travel  time  for  potential  water  contaminants  to  move  to  water 

24  resources?"    It  is  only  if  the  answer  to  this  question  is  yes  or  the  travel  time  is  short 

25  enough  to  be  of  concern  that  the  remaining  two  questions  are  important.    This  is  not 

26  to  suggest  that  the  movement  of  contaminants  in  water  out  of  the  active  rooting  zone 

27  is  not  of  concern.    It  is  simply  to  suggest  that  under  these  circumstances,  it  does  not 

28  constitute  a  risk  of  water  contamination.    This  accumulation  of  a  contaminant  under  dry 

29  conditions  over  the  long  term  resembles  the  condition  of  geologic  nitrate  described  by 

30  Follett  and  Walker  (1989). 

31  2.    what  kinds  of  potential  contaminants  are  present  and  in  what  amounts? 

32  At  a  detailed  level,  the  answer  to  this  question  is  quite  complex  and  includes  sources 

33  such    as   atmospheric    deposition,    release    from    forms    stored    in    soil    and    wetland 

34  components  as  well  as  the  incremental  inputs  used  in  the  course  of  land  husbandry.    In 
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1  this  discussion,    the   indicator   is  focused   on  the  specific   practices   associated    with 

2  agricultural  land  use.    It  may  be  thought  of  as  a  partial  budgeting  approach.    In  general 

3  then,  the  kinds  and  amounts  of  potential  contaminants  are  directly  related  to  the  kind 

4  of  crop  grown  and  the  type  of  agricultural  system.    From  the  standpoint  of  nutrients, 

5  there  are   three  general   sources;    namely,   commercial    fertilizer,   decomposing    crop 

6  residues  and  animal  manures.     The  amounts  can  be  inferred  from  the  current  crop 
production    recommendations,    the    amounts    harvested    in    crop    yields,    and    survey 

8  information  about  the  amounts  produced  or  sold.   For  pesticides,  the  kinds  and  amounts 

9  can  be  inferred  from  the  types  of  crops,  current  recommended  herbicide  usage  and 

1 0  surveys  of  actual  use. 

1 1  3.     what  is  the  quantity    of  water  in  excess  of  that  used  for  crop  production   and 

12  evaporation?    and,  what  is  the  fate  of  the  excess  water? 

13  These  are  probably  the  most  difficult  and  the  most  critical  questions  because  they 

14  control  the  extent  and  direction  of  movement  and  the  quantity  of  water  determines  the 

15  dilution  factor  for  all  potential  contaminants. 

16  3.2  ASSUMPTIONS 

1 7  At  all  but  the  lowest  levels  of  the  hierarchy,  reasonable  starting  assumptions  for  agricultural 

18  areas  include: 

19  1.      an  approximate  balance  between  nutrient  additions  and  removals  within  the  limits  of 

20  nutrient  use  efficiency. 

21  2.      an  indicator  which  estimates  the  impact  of  agricultural  (anthropogenic)  practices  can  be 

22  calculated  using  a  partial  budgeting  approach  which  considers  those  inputs  and  outputs 

23  directly  influenced  by  management  such  as  the  application  of  fertilizers,  manures  and 

24  pesticides,  and  the  removal  of  harvested  grain  and  stover. 

25  3.      trends  in  pesticide  use  can  be  estimated  from  i)  a  knowledge  of  the  crops  grown  and 

26  the  recommended  pesticides  associated  with  their  production,  and  ii)  the  general  level 

27  of  pesticide  usage  (e.g.,  county  level  surveys). 

28  4.      for  "average"   conditions  of  soil,  topography  and  drainage  it  is  possible  to  estimate 

29  IROWC  by  crop  (e.g.,  compared  to  non-crop  or  unmanaged  conditions). 

30  5.      at  lower  levels  of  the  hierarchy  it  is  possible  to  improve  the  IROWC  estimates  by 

31  incorporating   actual  conditions   of  soil  texture,  topography,   drainage   etc  and  more 

32  detailed  climatic  data. 
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1  6.      it  is  necessary  to  estimate  the  extent  to  which  "friendly"  practices  (e.g.,  conservation 

2  tillage,  buffer  strips)  reduce  the  IROWC  from  that  observed  under  "average"  conditions. 

3  3.3  A  PROPOSED  BUDGET  MODEL  FOR  IROWC 

4  This  section  describes  a  model  for  IROWC  based  on  a  budget  approach.    The  model  is  viewed 

5  as  a  component  of  a  broader  methodological  approach  to  IROWC,  which  will  integrate  the 

6  budget  approach  with  other  tools  (e.g.,  GIS)  and  approaches  (e.g.,  spatial  analysis  of  risk).  The 

7  essential  function  of  the  model  is  to  compare  the  flow-weighted  mean  concentration  of  each 

8  potential  contaminant  to  an  appropriate  standard  (e.g.  drinking  water).     The  risk  of  water 

9  contamination   from  non-crop  and  unmanaged   sources  is  acknowledged,    but  the  approach 

10  .        considers  only  the  added  risk  of  water  contamination  resulting  from  agricultural  activities. 

1 1  The  Potential  Contaminant  Concentration  (PCC)  consists  of  two  parts:  i)  determination  of 

12  .the  Potential  Contaminant  Present  (PCP)  in  mg/ha,  and  ii)  calculation  of  Excess  Water  (EW) 

13  in  1/ha,  representing  the  available  quantity  of  surplus  water.    Thus: 

14  PCC.  =  PCP  (mg/ha)  3.3    (1) 

15  EW  (1/ha) 

16  c  =  contaminant  type  (e.g.,  nitrate,  atrazine) 

17  PCP  =  f(climate,  crop,  input  rates,  soil,  management,  etc.) 

18  EW  =  f(precipitation,  evaporation,  crop  moisture  use  (yield),  soil). 

19  At  each  level  of  the  hierarchy,  the  differentiating  characteristics  and  associated  factors  are 

20  incorporated  into  the  determination  of  the  amount  or  potential  contaminant  present  (PCP).   For 

21  example: 

22  PCP5  =        f(crop  type,  crop  area,  nitrogen  in  fertilizer  and  manure,  nitrogen  in 

23  harvested  crop,  herbicide  class  (e.g.,  triazine),  annual  climate  normals). 

24  The  calculated   value   is  a  relative   indication  where  the  non-crop  or 

25  unmanaged  state  is  0  and  risk  increases  in  proportion  to  nitrogen  content 

26  of  the  harvested  crops,  and  level  of  pesticide  use. 

27  PCP3  =         f(all  information  and  constraints  from  higher  levels  plus  crop  rotations. 

28  soil    texture,    subsurface    drainage,    conservation    practices,    daily    and 

29  monthly  climate,  post-harvest  precipitation,  specific  pesticide,  etc.).   The 

30  calculated    value    is    the    concentration     in    mg/ha    for    the    specific 

3 1  contaminant. 

32  where  5  and  3  correspond  to  levels  of  the  hierarchy.    Excess  water  is  determined  by  estimating 

33  the  flows  in  the  hydrological  cycle  in  increasing  detail  at  lower  levels. 

11 
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1  IROWC  is  characterized  as  the  ratio  of  the  average  PCC  to  the  concentration  allowable  to 

2  maintain  a  desired  standard  (e.g.,  drinking  water).    Thus: 

3  IROWCc    =     PCC.JIlgZl 3.3    (2) 

4  Maximum  Allowable  concentration  (mg/l)c 

5  IROWCComposite  =     Enc  =  ,    IROWQ  3.3    (3) 

6  Ideally,  this  calculation  would  be  carried  out  for  each  (precipitation)  event  which  constitutes 

7  a  risk  for  water  contamination.    In  practice,  this  will  only  be  possible  at  the  most  detailed  level 

8  of  assessment.     The  hierarchial  levels  at  which  the  indicator  is  required  (see  Table  1  in  the 

9  concept  paper)  determine  the  kinds  of  data  available  and  the  detail  and  precision  which  can  be 

10  achieved. 

1 1  In  hierarchical  fashion,  we  can  start  at  the  highest  level  and  determine  the  constraints,  and  then 

12  consider  lower  levels  within  these.    For  example,  it  is  feasible  to  characterize  IROWC  at  the 

13  ecoregion  or  ecodistrict  level  (levels  4-5).     The  proportion  of  the  ecoregion  or  ecodistrict 

14  occupied  by  various  crops  can  be  determined.    At  these  levels,  only  a  relative  indicator  based 

15  on  annual  average  conditions  can  be  estimated.    It  is  sensitive  to  changes  in  crop  type,  crop 

16  area,   livestock   density  and  monthly  climate.      It  is  not  feasible   to  incorporate   "friendly" 

17  practices  because  it  is  not  possible  to  associate  them  with  specific  crops  (see  Table  1  in  the 

18  concept  paper). 

19  At  lower  levels  (levels  1-3),  the  indicator  will  more  nearly  reflect  actual  conditions.    It  will 

20  require  many  more  types  of  data  in  much  greater  detail  to  incorporate  spatial  variability  in  soil, 

21  topography  and  drainage;  the  degree  of  adoption  of  friendly  practices  and  more  specific  details 

22  of  climate.    Temporal  sensitivity  is  also  better  at  lower  levels  allowing  for  calculations  of  risk 

23  on  a  monthly  or  event  basis.    With  the  calculation  of  IROWC  on  the  basis  of  soil  units  (e.g. 

24  SLC  for  level  3,  detailed  maps  for  level  2  and  on-site  soil  characterization   for  level  1),  it 

25  should  also  be  possible  to  partition  the  risk  of  contamination  into  surface  or  subsurface  water 

26  and  possibly,  tile  drainage  water.    In  addition,  with  movement  from  highest  to  lowest  levels, 

27  the  probability  of  coincident  occurrence  (the  certainty  that  the  specific  crop  is  grown  on  a 

28  known  soil  type  using  particular  tillage  practices,  etc)  increases,  thereby  enhancing  the  level 

29  of  confidence  for  the  estimate. 

30  4.  EXAMPLE  APPLICATIONS  OF  THE  PROPOSED  METHODOLOGY 

31  The  following  sections  have  been  prepared  at  the  request  of  the  indicator  team  to  demonstrate 

32  how  the  proposed  methodology  would  characterize  the  risk  of  water  contamination.     A  full 

33  implementation  of  the  proposed  methodology  is  a  substantial  undertaking  requiring  consultation 

34  and  input  from  many  experts  from  across  Canada  and  internationally.    In  addition,  it  is  to  be 

12 
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1  expected  that,  as  details  of  the  methodology   are  worked  out,  there  will  be  gaps  in  our 

2  knowledge  and  understanding  which  require  original  research,  and  gaps  in  the  available  data 

3  which  will  require  alternative  data  sources  or  data  collection  activities.     This  could  include 

4  coupling  other  agri -environmental  indicators  with  IROWC,  or  direct  incorporation  of  their  data. 

5  Consequently,     THE    SAMPLE    IROWC     CALCULATIONS     SHOWN    HERE    ARE    A 

6  DEMONSTRATION    OF   THE   METHODOLOGY    ONLY.      THEY    HAVE   NOT   BEEN 

7  EXTENSIVELY  CHECKED  AND  SHOULD  NOT  BE  TREATED  AS  RELIABLE. 

8  These  examples  do  however  provide  a  demonstration  of  the  kinds  of  information  which  are 

9  required,  the  kinds  of  resolution  which  can  be  achieved  and,  very  clearly,  the  levels  of  effort 

1 0  and  detail  which  are  needed  to  calculate  something  approaching  ACTUAL  indicators  of  risk 

1 1  of  water  contamination. 

12  4.1  IROWC  AS  A  RELATIVE  MEASURE  -  BROAD  LEVEL  INTERPRETATION 

13  Ideally,  the  proposed  methodology  would  be  implemented  at  all  hierarchial  levels.    In  reality, 

14  it  is  only  possible  to  determine  all  parameters  required  at  the  lower  levels.   At  the  higher  levels 

15  only  a  relative  IROWC  rating  is  possible  based  on  the  principles  outlined;  namely,  an  estimate 

16  of  the  amount  of  potential  contaminant(s)  available,  the  amount  of  excess  water  and  the  amount 

17  which  could  be  accommodated  within  the  standards  of  concern. 

18  For  the  broad  level  assessments  (levels  4  and  5),  the  IROWC  is  based  on  current  land  use 

19  practices  which  give  an  indication  of  the  level  and  kinds  of  inputs  and  climate  records  which 

20  provide  a  rough  estimate  of  the  quantity  of  water  in  excess  of  crop  requirements.    These  results 

21  could  be  presented  at  the  level  of  ecozone,  ecoregion  or  ecodistrict.    The  difference  would  be 
in  the  degree  of  resolution  of  the  classes  and  the  extent  to  which  specific  land  management 

23  practices  (such  as  conservation  tillage,  buffer  strips,  etc)  could  be  considered.    In  addition,  for 

24  the  ecodistricts  it  would  be  possible  to  look  at  climate  patterns  for  portions  of  the  year  as  well 

25  as  the  degree  of  variation  between  years. 

For  this  demonstration  example,  the  results  are  presented  at  the  ecodistrict  (ED)  level.  This 
level  and  all  larger  spatial  units  are  so  large  that  the  land  use  and  management  practices  cannot 
be  associated  with  specific  soils,  topography  and  drainage  systems.  Consequently,  it  is  not 
possible  to  partition  the  IROWC  between  surface  and  ground  water.  Variations  in  specific  land 
management  practices,  or  variation  throughout  the  year  or  between  years,  can  be  incorporated 
31  but  are  not  illustrated  in  these  examples. 

The  moisture  conditions  were  estimated  from  climatic  normals  data  for  1951-80.    These  data 
were  extracted  from  the  Land  Potential  Data  Base  (Kirkwood  et  al.,  1983)  and  associated  with 
34  ecodistrict  (ED)  and  SLC  polygons  by  taking  the  largest  area  of  overlap. 


IROWC  -  Methodology  (DRAFT) 

1  In  keeping  with  the  first  question  posed  (section  3.1.4)  the  relative  annual  excess  water  (REW) 

2  or  deficit  was  calculated  from  the  annual  precipitation  (P)  and  the  potential  evaporation  (PE), 

3  using  1951-80  climatic  normals  data,  as: 

4  REW  =  Annual  (P  -  PE)  4.1    (1) 

5  The  estimate  of  REW  was  determined  for  the  prairie  region  and  southern  Ontario  (Figure  2). 

6  In  the  prairies,  all  but  the  fringe  areas  are  characterized  by  REW  less  than  zero.    While  there 

7  will  be  some  movement  of  water  to  surface  recharge  areas  at  specific  times  of  the  year  and 

8  also  some  movement  to  ground  water,  at  a  broad  level  this  movement  is  negligible.    This  is 

9  consistent  with  the  general  landscape  patterns  of  the  great  plains  where  most  of  the  drainage 

10  is  internalized  into  prairie  sloughs;  the  water  tables  are  quite  deep  with  little  connection  to 

1 1  surface  recharge  areas  and  there  are  few  external   drainage  channels.     This  assessment    is 

12  consistent  with  the  observations  of  de  Jong  and  Kachanoski  (1987).   For  further  demonstrations 

13  of  the  proposed  methodology,  Southern  Ontario  is  used  as  the  sample  region. 

14  It  is  recognized  that  this  approach  to  calculating  annual  excess  water  is  not  without  problems. 

1 5  Actual  evapotranspiration  (AE)  is  only  a  fraction  of  PE  during  most  of  the  day  in  regions  such 

16  as  the  prairies.    This  means  that  the  above  approach  may  underestimate  annual  excess  water. 

17  However,  the  approach  is  maintained  because  it  is  difficult  to  obtain  data  for  AE  at  broad 

18  scales.    A  related  problem  is  basing  spatial  screening  on  a  single  criterion  (i.e.,  REW).    In  the 

19  prairie  region,  it  is  possible  to  have  both  a  small  REW  and  a  large  potential  contaminant 

20  concentration   (PCC).     The  criterion  of  EW  is  maintained  over  concentration   because  the 

21  fundamental  matter  of  concern  is  water  transport  of  the  contaminant.     Again,  this  does  not 

22  imply  that  the  steady  accumulation  of  a  contaminant  under  dry  conditions  and  over  the  long 

23  term  (e.g.,  geologic  nitrate)  is  not  important,  but  rather  that  its  potential  for  transport  to  surface 

24  or  ground  water  is  significantly  reduced. 


I 
I 
I 
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25  Ideally,  PCC  should  incorporate  variation  in  nutrient  uptake  efficiency  and  time  of  uptake  for m 

26  each  crop,  as  well  as  the  temporal  variability  of  excess  water.    Thus: 

27  PCC  =  ffcrop  type,  nutrient  uptake  efficiency,  time  of  uptake)  4.1    (2) 

28  excess  water  f(time  period)  x  crop  area 

29  Data  constraints  and  the  comprehensiveness  of  the  analysis  resulted  in  a  simplified  version  of^ 

30  this  equation  for  the  purposes  of  this  demonstration.    This  analysis  used  mean  annual  values 

31  for  P  and  PE  (1951-1980  climatic  normals)  to  estimate  REW,  and  a  mean  nutrient  uptake! 

32  efficiency  of  60%  for  all  crops  over  the  entire  growing  season. 

33  Based  on  the  above  rationale,  the  following  series  of  equations  were  used  to  estimate  the  broadB 

34  level  IROWC  for  Southern  Ontario: 
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IROWC  -  Methodology  (DRAFT) 

1  RPCCN  =  Zed  Crop  Area^  CTOp   x  Nitrogen  Harvested by  CTOpj  4.1    (3) 

2  REWx  total  crop  area 

3  RPCCM  =  EFP  Nitrogen  in  Manure  x  0.1  4.1    (4) 

4  REW  x  total  crop  area 

5  RPCCy,  =  ZFP  Crop  Area.,.,  frop    x  Average  rate  for  herbicides by  „„r  4.1    (5) 

6  REW  x  total  crop  area 

7  RPCCHa  -  RPCCH  x  fP-PE)M...  +  (P-PE),.„.  x  0.5  +  (P-PE)„„,  x  0.25  4.1    (6) 

8  (Annual  Precipitation  -  PE  from  1951-80  climatic  normals) 

9  RPCC  =  relative  potential  contaminant  concentration  for  nitrogen  in  the  harvested  crop  (N), 

10  long-term  nitrogen  from  manure  (M),  herbicide  (H),  and  herbicide  concentration  adjusted  for 

1 1  the  susceptible  time  (month)  for  transport  (Ha).  The  herbicides  atrazine  and  simazine  are  used 

12  for  demonstration  purposes. 

13  IROWC(level  4.5)  =  (RPCCN  +  RPCCM)/10  +  RPCCHa/0.01  4.1    (7) 

14  The  IROWC  provides  a  measure  of  the  relative  amount  by  which  the  current  standards  are 

15  exceeded.    The  current  drinking  water  standard  for  nitrogen  is  10  mg/1  and  the  standard  for 

16  atrazine  is  0.005  mg/1,  and  for  simazine  it  is  0.01  mg/1.   The  latter  standard  is  used  to  calculate 

17  IROWC  for  herbicides. 

18  Three  additional  considerations  were  included  in  the  these  equations: 

19  1.     Data  on  nitrogen  use  efficiency  plotted  against  proportion  of  optimum  yield  (Sander  et 

20  al.   1994)  shows  an  efficiency  of  60-40%  for  corn.     Data  from  Europe  (Juergens- 

21  Gschwind  1989)  also  suggests  a  nitrogen  use  efficiency  of  60%  and  further  indicates 

22  the  other  sink  terms  for  the  additional   nitrogen  (e.g.  denitrification   to  GHG,   soil 

23  biomass,  transport  off-site  by  water).    This  implies  that  the  quantity  available  for  other 


24  '  The  calculation  determines  the  average  quantity  of  nitrogen  harvested  per  hectare  of 

25  cropped  land  in  an  ecodistrict  based  on  the  following  levels  of  nitrogen  in  the  harvested  crop: 

26  Corn,  grain  and  silage  -  175  kg/ha;  soybeans  and  white  beans  -  150  kg/ha;  barley  -  45  kg/ha; 

27  spring  wheat  -  55  kg/ha;  winter  wheat  -  85  kg/ha;  fall  and  spring  rye  and  oats  -  35  kg/ha; 

28  alfalfa  and  other  tame  hay  -  200  kg/ha. 
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IROWC  -  Methodology  (DRAFT) 

1  pathways  is  approximately  equivalent  to  the  Nitrogen  Harvested  (equation  3).    This 

2  efficiency  is  assumed  for  all  crops. 

3  2.     OMAF  Publication  296  indicates  the  long-term  value  of  manure.   Generally,  50%  of  the 

4  total  manure  nitrogen  (75%  for  poultry)  applied  in  the  spring  is  available  in  the  year  of 

5  application.    In  the  second  year  approximately  10%  of  the  remaining  nitrogen  becomes 

6  available,  5%  in  the  third  year.    For  this  calculation,  the  nitrogen  during  the  year  of 
application  is  included  in  the  estimate  based  on  harvested  crop.    For  subsequent  years 

8  a  factor  of  1 0%  is  used  to  account  for  the  long-term  value  and  also  for  the  nitrogen  lost 

9  from  manure  applied  (equation  4). 

10  3.     For  herbicides,  equation  6  is  adjusted  by  the  proportion  of  the  excess  moisture  in  May, 

1 1  June  and  July  as  these  are  the  times  the  herbicide  will  be  most  susceptible  to  transport. 

12  These  monthly  figures  are  further  weighed  (1.0,  0.5  and  0.25  respectively)  to  account 

13  for  declining  susceptibility  to  movement  off-site.    After  this  time  it  is  assumed  to  have 

14  been  degraded  by  natural  processes.    Atrazine  and  simazine  are  aggregated  as  triazines, 

15  and  specific  properties  of  individual  herbicides  (e.g.,  solubility,  half-life,  partitioning 

16  coefficient)  are  not  considered  at  levels  4-5. 

17  Key  features  demonstrated  in  the  example  application  at  the  level  4-5  assessment  are  trends  and 

18  one  option  for  calculating  a  composite  indicator. 

19  4.1.1  Temporal  variation 

20  Figure    3    shows    the    change    between    1991    and    1981    in   relative    potential    contaminant 

21  concentration   for  nitrogen  (based  on  levels  harvested  in  the  crop)  for  crop  distribution  as 

22  reported  in  the  census  of  agriculture.    These  results  are  reported  in  units  of  mg/1  or  ppm  with 

23  the  range  going  from  a  decrease  of  5  or  less  to  an  increase  of  5  or  more.    In  general  these 

24  results  are  quite  encouraging  for  much  of  the  province  showing  no  change  or  a  slight  decline. 

25  The  greatest  increase  is  indicated  in  the  ecodistrict  which  occupies  north  Middlesex  and  the 

26  western  portion  of  Huron  Counties.     In  addition,  there  is  some  indication  of  an  increase  in 

27  Eastern  Ontario  and  Essex.  Lambton  counties  and  the  fringe  area.   The  map  illustrates  that  this 

28  kind  of  calculation   can  indicate  some  changes  which  are  potentially  of  the  same  order  of 

29  magnitude  as  the  drinking  water  standard. 

30  4.1.2  Composite  indicator 

31  The    ratio    of  the    relative    potential    contaminant    concentration    (RPCC)    to   the   maximum 

32  concentration  allowed  by  current  standards  provides  an  indication  of  risk.    If  all  water  quality 

33  standards  have  the  same  level  of  importance  (or  risk  if  exceeded)  then  it  is  possible  to  provide 

34  a  composite  indication  of  risk  by  combining  the  risks  for  all  separate  potential  contaminant 

35  components.    This  aspect  is  described  by  equation  7  and  illustrated  for  Figures  4a,  4b  and  4c. 


IROWC  -  Methodology  (DRAFT) 

1  Figure  4a  shows  the  relative  IROWC  associated  with  nitrogen  from  crop  production  and  also 

2  the  residual   from  previous   manure   applications.      Figure  4b  shows   the  relative   IROWC 

3  associated  with  triazine  applications.    Figure  4c  shows  the  results  of  the  combination  of  these 

4  component  IROWC  estimates. 

5  4.2  IROWC  AS  AN  ACTUAL  MEASURE  -  ASSESSMENTS  AT  LOWER  LEVELS  (3-1) 

6  As  we  move  down  the  hierarchy,  the  resolution  and  accuracy  of  the  IROWC  should  improve. 
.  7  In  addition   to  the  dominant   factors   affecting   contamination   at  each   higher   level    in  the 

8  hierarchy,  an  assessment  of  risk  at  lower  levels  incorporates  greater  detail  through  refinement 

9  of  parameters  such  as  climate    (e.g.,  monthly  values  of  P-PE  (or  P-AE)  to  identify  months 

10  when  off-site  movement  is  most  likely),  nutrient  inputs  (e.g.,  nitrogen  uptake  pattern  of  specific 

1 1  crops  over  the  growing  season  to  determine  when  surplus  nitrogen  is  most  likely  available). 

12  tillage  practices  (e.g.,  conventional,  conservation,  no-till)  and  their  effect  on  excess  water  and 

1 3  surplus  nutrients,  and  processes  such  as  subsurface  tile  drainage  to  further  partition  surface  and 

14  ground  water  components. 

15  An  assessment  of  IROWC  at  lower  levels  incorporates  these  and  other  detailed  parameters  in 

16  order  to  estimate  an  actual  concentration  of  a  contaminant  in  water.    In  most  cases,  this  will 

17  result  in  requirements  for  additional  data  and  lead  to  increasingly  complex  analyses.    It  was 

18  not  possible  to  develop  a  complete  example  demonstration  of  IROWC  at  lower  levels.    The 

19  following  sections  provide  OUTLINES  of  how  COMPONENTS  of  IROWC  are  assessed  at 

20  lower  levels.    They  also  DEMONSTRATE  the  kinds  of  data  which  must  be  synthesized.    The 

21  results  are  closer  to  ACTUAL  BUT  THEY  ARE  STILL  RELATIVE. 

22  Several   attributes   of  IROWC   are   selected   to  demonstrate   the  type  of  parameters   to  be 

23  considered  at  lower  levels.    The  geographic  focus  is  on  watersheds  in  south-western  Ontario 

24  where  flows  and  sediment  yield  have  been  monitored  and,  one  small  watershed  in  particular, 

25  Bamberg  Creek. 

26  4.2.1  Bamberg  Creek  Watershed 

27  This  watershed  is  a  catchment  located  in  Wilmot  and  Wellesley  Townships  in  the  Regional 

28  Municipality  of  Waterloo.    Agriculture  occupies  about  81%  of  the  watershed's  land  area  in  a 

29  mix  of  farm  types  (e.g.,  cash  crop,  livestock,  mixed  farming).   A  recent  study  of  the  watershed 

30  was  undertaken  to  assess  the  state  of  the  agricultural  land  base  relative  to  its  contributions  to 

31  non-point    source    pollution    (Ecologistics    Limited    1993).      The   study    integrated    available 

32  information  on  soils,  land  use  and  farm  inputs  (fertilizers,  manures  and  pesticides)  into  a  GIS 

33  framework  to  produce  maps  (1:50,000)  of  nutrient  and  pesticide  loading,  and  ground  water 

34  contamination  potential.    The  approach  used  approximates  the  preliminary  steps  of  a  budget 

35  approach.    Furthermore,  data  on  soils,  land  use,  and  inputs  of  nitrogen  and  pesticides  used  in 

36  this  study  indicate  the  information  base  available  for  a  level  3  assessment. 
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IROWC  -  MethodologN   (DRAFT) 

1  Nitrogen  inputs  were  estimated  using  fertilizer  nitrogen  recommendations  for  the  various  field 

2  crops  grown.     These  were  converted   to  a  unit  area  loading   for  each   land  use  category 

3  (assuming  different  proportions  of  specific  crops  in  various  land  use  systems)  and  weighted  by 

4  the  proportion  of  each  land  use  assumed  to  have  received  fertilizer.   Estimates  of  nitrogen  from 

5  manure  were  based  on  animal  units  (census  data),  manure  production  per  animal  unit,  and 

6  average  nitrogen  concentration  of  manure.     These  values  were  adjusted  to  account  for  an 

7  estimated  nitrogen  loss  of  50%  during  storage,  handling  and  application.    Areal  proportions  of 

8  manure  application  per  land  use  system  were  assumed  to  determine  unit  area  loadings.    The 

9  nitrogen  values  for  fertilizer  and  manure,  and  total  nitrogen,  are  shown  by  land  use  in  Table 

10  2.    These  data  were  used  to  produce  a  map  of  nitrogen  loading  (Figure  5). 

11  Using  1986  census  information,  and  data  from  the  1988  Ontario  pesticide  survey,  the  same 

12  basic  approach  was  used  to  estimate  pesticide  loadings  by  land  use  category.    Pesticides  were 

13  not  differentiated  by  type  in  this  study.    Crop  groups  were  defined  (row  crops,  small  grains. 

14  hay  and  improved  pasture)  and  assumptions  made  to  estimate  unit  area  pesticide  loading  for 

15  each  group  and  distributed  over  the  land  use  categories.    Calculated  pesticide  loadings  by  land 

16  use  category  are  shown  in  Table  2. 

17  Finally,  a  map  of  soil  transmissivity  (i.e.,  ability  of  surficial  soil  (<lm)  to  transmit  pollutants) 

1 8  was  superimposed  onto  maps  of  nitrogen  and  pesticide  loading  to  display  the  relative  potential 

19  for  ground  water  contamination  from  these  contaminants. 

20  4.2.2  Extending  the  Research 

21  The   focus  of  the  Bamberg   Creek   study   was  on  characterizing   the  amounts   of  potential 

22  contaminants  present.     From  the  standpoint  of  an  IROWC,   there  are  additional   data  and 

23  information  requirements.    The  analysis  presented  results  on  an  annual  basis  with  no  attempt 

24  to  deal  with  changing  conditions  throughout  the  year.     There  was  no  attempt  to  develop 

25  complete  budgets  for  nutrients  or  pesticides,  and  the  complex  and  difficult  task  of  partitioning 

26  water  between  surface,  tile  flow  and  ground  water  was  not  addressed.  There  are  no  easy  ways 

27  of  dealing  with  these  concerns  but  the  following  sections  indicate  some  possible  approaches 

28  and  data  sources  which  can  be  used. 

29  1.  Temporal  Distribution  of  Surplus  Nitrogen  and  Excess  Water.    Nitrogen  is  present  in 

30  the  soil  throughout  the  year  and  is  normally  added  as  fertilizer  early  in  the  growing  season. 

31  The  uptake  of  nitrogen  by  plants  increases  over  the  growing  season.    The  pattern  of  nitrogen 

32  uptake  can  be  analyzed  in  relation  to  surplus  moisture  distribution  to  determine  the  amount  of 

33  mobile  nitrogen  in  the  soil  profile. 

34  Ontario  data  on  plant  uptake  of  nitrogen   is  available   for  some  crops  (M.  Miller,   personal 

35  communication).    The  quantity  of  nitrogen  taken  up  by  corn  (grain  and  stover)  is  illustrated  in 

36  Figure  6.    Uptake  commences  after  the  first  month  of  growth,  increasing  steadily  between  32- 
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Table  2  Nitrogen  and  Pesticide  Loading  by  Land  Use  in  Bamberg  Creek  Watershed 


Land  use 

Fertilizer  N     Manure  N 

Total  N 

Pesticides 

-  kg  N  ha"'  - 

(kg  a.i.  ha"1) 

Continuous  crop 

226                   92 

318 

1.60 

Corn  system 

190                  133 

323 

1.08 

Mixed  system 

165                   162 

327 

0.73 

Grain  system 

102                   39 

141 

0.40 

Hay  system 

126                  342 

468 

0.33 

Pasture  system 

54                   235 

289 

0.10 

Grazing  system 

0                    169 

169 

source:  data  from  Ecologistics  Limited  (1993) 

1  112  days  after  planting  (DAP),  and  levelling  off  somewhat  following  silking  (112  DAP).   This 

2  pattern  suggests  that  all  of  the  applied  nitrogen  is  surplus  during  the  early  stage  of  plant 

3  growth.    During  subsequent  growth  stages,  nitrogen  uptake  increases  steadily  (i.e.,  the  amount 

4  of  surplus  nitrogen  decreases).     Over  the  entire  growing  season,  total  nitrogen  uptake  is 

5  approximately  200  kg  ha"1. 

6  The  temporal  distribution  of  nitrogen  uptake  can  be  compared  to  the  seasonal  or  monthly 

7  distribution  of  surplus  moisture  in  order  to  estimate  the  quantity  of  nitrogen  at  risk  of  leaching. 

8  An  actual  estimate  can  be  derived  using  data  on  monthly  or  daily  climatic  normals  (e.g., 

9  precipitation,  potential  evapotranspiration)  for  a  representative  climate  station.  An  example  for 

10  the  Bamberg  Creek  watershed  using  climate  data  from  Elora  Research  Station  is  shown  in 

1 1  Figure  6.    There  is  a  mean  annual  water  surplus  of  367  mm,  but  the  distribution  is  such  that 

12  a  water  deficit  exists  during  the  summer  months  (June  -  September).    This  suggests  little  risk 

13  of  nitrogen  movement  during  the  growing  season.    Even  30  days  after  planting,  when  nitrogen 

14  uptake  is  nil  (i.e.,  surplus  nitrogen  is  high),  the  risk  is  low  because  the  onset  of  water  deficit 

15  coincides  with  planting  of  corn  and  fertilizer  application  (i.e.,  mid-late  May).     Of  course 

16  individual  rainfall  events  may  contribute  to  movement  of  surplus  nitrogen  throughout  the 

17  growing  season,  as  well  as  residual  nitrogen  in  the  soil  at  post-harvest.     This  would  require 

18  further  detailed  analyses  using  daily  climate  daily. 

19  Crop  uptake    is  an  important,    but  incomplete,    measure    of  recovery    of  applied   nitrogen. 

20  Nitrogen  is  also  partitioned  in  crop  residue  (surface  stubble,  roots)  and  in  the  soil  as  residual 

21  nitrogen.    Like  crop  uptake,  these  also  vary  by  crop  type  and  rate  of  nitrogen  application.    It 

22  will  be  necessary  to  incorporate  the  changes  in  nutrients  and  water  throughout  the  year  to 

23  produce  a  level  3  IROWC. 
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1  2.  The  use  of  Computer  Simulation  Models  to  Characterize  the  Dynamics  of  specific 

2  Pesticides  in  Soil  Water. 

3  Another  way  to  incorporate  the  approach  and  results  of  Bamberg  Creek  watershed  study  into 

4  an  assessment  of  IROWC  at  level  3  is  to  determine  the  fate  of  specific  pesticides  using  models 
which  estimate  residual  pesticide  movement  out  of  the  root  zone.    The  study  aggregates  data 

6  on  pesticide  use  but  does  not  differentiate  pesticides  by  type  or  rate  of  application.    Pesticides 
vary  in  their  physical  and  chemical  properties,  and  in  their  soil-water  interactions.    An  analysis 

8  which  considers  various  chemical  transformations  and  partitioning  of  pesticide  residue  into 

9  differing  sinks  is  best  approached  by  use  of  a  pesticide  fate  model  (e.g.,  LEACHM,  GLEAMS, 

10  PRZM).    The  spatial  information  on  pesticide  loading  in  Bamberg  Creek  watershed  could  be 

1 1  used  to  target  specific  sites  (field,  farm)  for  more  detailed  risk  analyses  of  pesticide  movement. 

12  In  addition  to  site-specific  applications,  it  is  also  possible  to  apply  pesticide  fate  models  at  the 

13  watershed  scale.   A  recent  study  by  de  Jong  et  al.  (1994)  developed  and  applied  a  methodology 

14  to  investigate  migration  of  atrazine  through  the  soil  profiles  (<90  cm)  of  the  Grand  River 

15  watershed.    This  watershed  contains  the  Bamberg  Creek  catchment.    The  study  developed  an 

16  integrated  methodological  framework  consisting  of  1)  modified  submodels  of  LEACHM  to 

17  describe   soil  water  flow  (LEACHW)   and  atrazine  movement  (LEACHP),   2)  pedotransfer 

18  functions  based  on  soil  texture  and  organic  carbon  content,  3)  geostatistical  analyses  to  extend 

19  soil  and  climate  data,  and  model  output,  from  point  to  areal  format,  and  4)  GIS  to  produce 

20  maps  of  atrazine  loading  in  the  Grand  River  watershed  (1:1,000,000). 

2 1  Results  of  this  study  include  model  predictions  of  annual  atrazine  loading  over  time  for  several 

22  soils  (Figure  7).    Loadings  vary  by  soil  type,  but  tend  to  stabilize  over  time  (after  5-8  years) 

23  for  all  soils.    Two  of  the  soil  types  are  found  in  the  Bamberg  Creek  watershed  (Huron  clay, 

24  Fox  sand).    Atrazine  movement  in  the  Huron  soil  is  characterized  by  an  incremental  increase 

25  in  loading,  stabilizing  after  year  three.     Presumably,  atrazine  applied  to  this  soil  is  readily 

26  degraded  or  its  movement  out  of  the  soil  profile  is  impeded  (de  Jong  1994).    The  Fox  sand  is 

27  distinguished  by  a  relatively  high  atrazine  loading,  posing  a  greater  risk  of  water  contamination. 

28  Seasonal  distribution  of  atrazine  loading  is  not  known  since  model  simulations  were  undertaken 

29  on  an  annual  basis.    Other  factors  (e.g.,  information  on  land  use,  crop  rotations,  and  changes 

30  in  water  table  depth)  were  excluded  in  the  analysis,  but  the  methodology  has  considerable 

31  potential  for  analyzing  the  risk  of  pesticide  contamination  of  ground  water,  and  for  projecting 

32  the  impact  of  changes  in  land  use  and  pesticide  management  at  the  watershed  scale. 

33  3.  Partition  Excess  Water  into  Surface  Tile  Drainage  and  Ground  Water 

34  A  final  suggestion  for  building  on  the  Bamberg  Creek  watershed  study  to  calculate  a  level  3 

35  IROWC  assessment  is  to  partition  the  excess  water  component  into  surface  or  ground  water. 

36  Ultimately,  detailed  analyses  to  quantify  hydrologic  pathways  (e.g.,  runoff,  base  flow,  deep 
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1  percolation)  is  best  undertaken  by  hydrology  models  and  use  of  long  term  data  on  climate, 

2  stream  flow,  etc.    However,  it  may  be  possible  to  use  other  data  to  approximate  the  water 

3  partitioning.    These  alternative  data  sources  include: 

4  -     estimates    based   on   stream   flow  to  quantify   amount   of  water   available   to   move 

5  contaminants, 

6  -     analysis  of  land  resource  and  management  factors  which  influence  moisture  partitioning, 

7  and 

8  -     output  from  the  hydrology  submodel  of  models  already  used  to  study  sediment,  nutrient 

9  or  pesticide  movement. 

10  The  Monitoring  and  Systems  Branch,  Ontario  Region,  Environment  Canada,  collects  data  on 

1 1  stream   flow  and  sediment   levels   for  a  number  of  watersheds    in  south-western    Ontario, 

12  including  the  watershed  where  Bamberg  Creek  is  located.    If  we  assume  that  the  ground  water 

13  system    is  approaching    equilibrium    then   stream    flow  will   approximate    the   excess   water 

14  available.    This  measurement  does  not  provide  information  on  the  portion  of  the  water  which 

15  is  surface  runoff  or  the  portion  which  has  moved  through  the  soil.   Figure  8  shows  the  location 

16  of  the  9  watersheds  where  stream  flow  data  were  available.   Table  3  compares  the  excess  water 

17  as  determined  from  climatic  normals  data  to  the  1972-1991  mean  annual  stream  flow.   Values 

18  of  the  stream  flow  estimate  suggest  that  the  actual  quantity  of  excess  water  is  larger  and  more 

19  variable  than  the  estimate  based  on  climate.    These  results  are  to  be  expected  because: 

20  -  the  climate  data  has  been  generalized  over  relatively  large  areas,  and 

21  -  actual  evapotranspiration  is  always  less  than  potential. 

22  The  partitioning  of  water  between  the  soil  surface  and  subsurface  is  determined  by  many 

23  factors.     Rigorous  estimates  are  quite  difficult  to  make.     There  are,  however,  a  number  of 

24  factors  of  soil,  landscape  and  management  which  can  be  used  to  provide  a  more  qualitative 

25  assessment.    The  soil  polygons  shown  on  Figure  8  partition  the  watersheds  into  different  areas 

26  depending  on  the  land  resource  properties.   Characteristics  such  as  texture,  slope  and  landform 

27  have  a  strong  influence  on  water  pathways.   Sandy  textured  soils  are  more  permeable  and  allow 

28  more  flow  through  the  soil  and  less  runoff.     Areas  of  greater  slope  steepness  and  more 

29  dissected    landform    are   conducive    to   higher    proportions    of  runoff.      These    factors    are 

30  summarized  in  Table  4  for  those  soil  polygons  which  are  more  than  50%  within  watershed 

3 1  boundaries. 

32  Other  factors  which  influence  the  water  partitioning  are  the  proximity  to  natural  or  artificial 

33  watercourses.    The  stream  network  shown  in  Figure  8  was  obtained  from  Environment  Canada. 

34  An  analysis  was  done  to  determine  the  proportion  of  land  adjacent  to  streams.   Land  which  was 

35  within  200  meters  of  a  stream  was  considered  to  be  in  close  proximity.    The  rationale  for  the 

36  200  meter  buffer  was  that  field  boundaries  tend  to  impede  surface  runoff  and  the  traditional 

37  dimensions  for  fields  in  Ontario  are  40  rods  square  or  200  x  200  meters.   The  stream  proximity 

38  was  calculated  for  each  soil  landscape  polygon  and,  as  shown  in  Table  4,  it  varied  substantially 

39  across  the  watershed  area. 
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Table  3  Comparison  of  Annual  P-PE  and  Annual  Stream  Flow  in  Selected  Watersheds 


Watershed     Station  No.  P-PE  (mm)  steam  flow  (mm) 


1 

02GB001 

220  -  320 

384 

2 

02GE003 

220  -  250 

422 

3 

02GG007 

160  -  320 

323 

4 

02GC002 

220  -  230 

362 

5 

02GC018 

220  -  230 

391 

6 

02GC026 

220  -  230 

449 

7 

02GC021 

220  -  230 

512 

8 

02GC007 

220  -  230 

398 

9 

02GB007 

220  -  230 

312 

Data  Sources: 

DP-PE  based  on  19. 

51-80  climatic  norma 

Is  (Kirkwo* 

2)Runoff  based  on  1972-91  observation  ofhydrologic  stations,  Monitoring  and  System  Branch, 
Environment  Canada. 


1  One  of  the  more  difficult  aspects  of  water  partitioning  and,  consequently,  the  destination  of 

2  potential  contaminants  is  flow  through  artificial  subsurface  drains.   Water  moving  through  this 

3  pathway  will  leach  materials  from  within  the  soil  profile  but  then  is  shunted  to  surface  waters 

4  rather  than  contributing  to  base  flow  or  ground  water  recharge.    While  the  soil  taxonomy  and 

5  natural  drainage  state  provide  some  basic  indications  of  presence  or  absence  of  tile  drains, 

6  maps  of  actual  tile  installations  must  be  used  to  provide  an  realistic  estimate  of  the  proportion 

7  of  drained  land.    Maps  of  tile  drainage  locations  were  available  for  portions  of  the  watersheds 

8  in  Oxford,  Middlesex  and  Elgin  counties.     Table  4  shows  these  estimates  where  data  were 

9  available  and  illustrates  the  substantial  degree  of  variation  to  be  expected. 

10  These  factors  illustrate  the  complexity  of  factors  associated  with  the  task  of  partitioning  water 

1 1  flows  in  the  environment.    Generally,  models  will  be  required  to  provide  a  realistic  integration. 

12  These  models  may  be  primarily  for  the  hydrologic  cycle  or  they  may  model  nutrient  and 

13  pesticide  flows  in  addition  to  water.    For  example,  some  hydrologic  data  should  be  available 

14  for  the  Bamberg  Creek  watershed  because  the  CREAMS  model  generally  partitions  runoff  and 

15  infiltration  (using  the  SCS  runoff  curve  number),  and  uses  a  water  balance  routing  to  account 

16  for  storage  in  soil  water,  uptake  by  plants,  and  deep  drainage. 
17 

18  Accurately  partitioning  excess  water  into  surface  and  ground  water  components  is  a  challenging 

19  hydrological  task  and  an  important  area  of  further  work. 
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Table  4     Some  Land  Resource  and  Management  Factors  which  Affect  the  Partitioning  of 
Excess  Water  for  Selected  Watersheds  of  Southern  Ontario  (data  compiled  at  SLC  level) 


Water  SLC 

Surface 

Slope 

Local 

Stream 

Prop. 

of 

shed 

poly 

#    Texture 

(%) 

Landform 

Proximity* 

Land  tiled 

1 

0078 

Silt  loam 

1  -  3 

Level 

10.6 

n/a 

0073 

Clay  loam 

1  -  3 

Undulating 

8.0 

n/a 

0452 

Clay  loam 

1           *"> 

1  -  j 

Undulating 

1.3 

n/a 

0066 

Clay  loam 

i      -> 
l  -  j 

Level 

5.6 

n/a 

0153 

Loam 

1  -  3 

Undulating 

4.8 

n/a 

0077 

Loam 

4  -  9 

Undulating 

3.6 

n/a 

0036 

Loam 

4  -  9 

Rolling 

7.4 

15.8 

0069 

Loam 

4  -  9 

Undulating 

6.6 

n/a 

0070 

Loam 

10  -  15 

Hummocky 

5.4 

n/a 

0471 

Loamy  sand 

1  -  3 

Undulating 

8.5 

n/a 

0470 

Loamy  sand 

1  -  3 

Undulating 

7.7 

n/a 

0049 

Loamy  sand 

1  -  3 

Undulating 

12.2 

4.5 

0067 

Sandy  loam 

1  -  3 

Undulating 

5.1 

n/a 

0068 

Sandy  loam 

4  -  9 

Rolling 

5.3 

n/a 

2 

0035 

Silt  loam 

4  -  9 

Undulating 

5.0 

47.7 

0028 

Silt  loam 

4  -  9 

Rolling 

12.0 

51.2 

0013 

Silt  loam 

1  -  3 

Undulating 

9.0 

26.2 

0038 

Clay  loam 

1  -  3 

Undulating 

8.9 

36.0 

0029 

Loam 

1  -  3 

Undulating 

7.9 

37.3 

0034 

Loam 

4  -  9 

Rolling 

6.3 

40.2 

3 

0026 

Silt  loam 

1  -  3 

Undulating 

8.7 

38.1 

0024 

Clay  loam 

1  -  3 

Level 

3.8 

44.9 

0020 

Sand 

1  -  3 

Undulating 

6.8 

4.9 

0025 

Sand 

1  -  3 

Undulating 

5.7 

19.2 

6 

0424 

Silty  clay  loam 

4  -  9 

Undulating 

13.3 

17.7 

0423 

Fine  sand 

1  -  3 

Undulating 

6.6 

8.3 

7 

0422 

Loamy  fine  sand 

1  -  3 

Undulating 

11.4 

n/a 

9 

0425 

Clay  loam 

4  -  9 

Ridged 

0.0 

15.5 

0145 

Loam 

1  -  3 

Undulating 

6.4 

n/a 

0031 

Loamy  sand 

1  -  3 

Undulating 

7.7 

8.5 

*  proportion 

of  land  (%)  within 

200  meters 

(40  rods)  of  streams  by  SLC  polygon 

ZJ> 
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1  4.2.3  IROWC  at  the  level  of  the  SLC  polygon 

The  discussions  of  data  for  a  level  3  calculation  of  IROWC  have  dealt  with  watersheds,  SLC 

3  polygons  as  shown  in  Figure  8  and  small  catchment  areas.    As  the  size  of  the  spatial  unit  gets 

4  smaller,  it  becomes  increasingly  difficult  to  assemble  comprehensive  consistent  data  over  large 

5  areas.    For  example,  it  was  possible  to  compile  census  of  agriculture  data  on  crop  distribution 

6  for  SLC  polygons.    Other  data  such  as  climate,  pesticide  usage  and  farming  systems  cannot  be 
directly  related  to  SLC  polygons.    The  climate  station  network  is  too  sparse  to  be  interpolated 

8  with  precision.    Pesticide  usage  is  recorded  on  a  substantially  larger  units  (counties)  in  classes. 

9  Farming  systems  can  be  compiled  from  census  but  requirements  for  confidentiality  preclude 

10  this  compilation  at  the  SLC  level. 

1 1  Based  on  the  crop  distribution,  the  relative  potential  contaminant  concentration  for  nitrogen 

12  based  on  crop  uptake  was  calculated  for  the  watershed  area.   Figure  9  shows  the  results  of  this 

13  calculation  for  1981  and  1991.   This  figure  shows  a  substantial  degree  of  change  over  the  time 

14  period  and  also  shows  the  increasing  complexity  which  results  as  the  spatial  units  become 

15  smaller  (i.e.,  approaching  an  actual  IROWC). 

16  5.  POTENTIAL  PILOT  AREAS  FOR  DETAILED  IMPLEMENTATION    OF  IROWC 

17  The  proposed   methodology   for  IROWC   must  be  implemented   and  tested   using  accepted 

18  scientific  principles  and  procedures.    As  it  is  implemented  it  will  be  important  to  review  and 

19  refine  the  data  and  scientific  basis  to  ensure  the  best  result.    Before  IROWC  can  be  accepted 

20  it  must  be  validated.     Testing  must  continue  on  an  ongoing  basis  to  ensure  accuracy  and 

21  relevance. 

22  At  the  general  level  (4-5),  IROWC  results  will  be  validated  by  comparison  to  data  from  the 

23  limited  monitoring  and  experimental  sites  across  Canada.     As  the  IROWC  methodology  is 

24  developed,  refined  and  finalized  it  will  be  appropriate  to  select  one  research  group,  with 

25  support  from  a  team  of  technical  experts,  to  assemble  the  data  and  implement  the  level  4-5 

26  IROWC.      Furthermore,   the  IROWC   level   4-5  results   will   be  compared   to  the   IROWC 

27  calculated  at  more  detailed  watershed  and  farm  levels. 

28  The  lower  level  (2-3)  IROWC  sites  become  important  locations  for  development   of  more 

29  accurate  indicators  and  also  as  comparison  sites  for  the  general  level  4-5  IROWC.     Their 

30  selection  should  be  based  on  a  variety  of  factors  including: 

31  -     locations  representative  of  the  soils,  land  use  and  management  and  geographical  areas 

32  of  agriculture  in  Canada 

33  -     sites  where  data  is  available  and  collected  on  an  ongoing  basis  to  calculate  the  moisture 

34  balance  and  the  budget  of  nutrients  and  chemicals  in  the  environment 

35  -     sites   where   past   and   ongoing   monitoring   activities   are   measuring   the   parameters 

36  necessary  to  test  moisture  and  budgeting  models. 
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1  -     in  light  of  current  resource  constraints,  it  will  likely  be  necessary  to  incorporate  the 

2  detailed  implementation  of  IROWC  into  closely  related  research  projects  which  are 

3  already  in  progress. 

4  The  following  is  a  list  of  some  potential  areas  for  detailed  IROWC  implementation: 

5  1.     the  Lower  Fraser  River  valley  of  B.C.  where  a  major  initiative  is  underway  to  estimate 

6  detailed  nutrient   budgets   for  specific   agricultural    locations.      A  focal   point   is  the 

7  Abbotsford  aquifer  where  a  substantial  project  has  been  carried  out  to  model  ground 

8  water  risks  from  intensive  livestock-based  agriculture. 

9  2.     the  Lethbridge  research  station  where  projects  relate  to  both  dryland  and  irrigated 

10  agriculture  and  monitoring  of  nutrient  and  pesticide  dynamics. 

11  3.     in  Ontario,  there  are  three  small  paired  watersheds  (Essex,  Kettle  and  Kintore)  where 

12  data  collection  and  monitoring  activities  have  been  carried  out  for  the  past  10  years. 

13  At  the  Essex  site  and  the  adjacent  Woodslee  Research  substation  detailed  hydrological 

14  measurements  are  being  taken.    At  the  Kintore  site,  there  is  a  large  project  currently 

15  under  way  to  characterize  water  movement  and  associated  transport  of  chemicals. 

16  4.     in  Eastern  Canada  there  are  a  series  of  small  watersheds  which  have  been  extensively 

17  characterized  and  monitored  in  association  with  the  International  Hydrological  decade. 

18  Since  characterization  of  the  hydrological  cycle  is  critical  to  the  success  of  the  IROWC 

19  methodology,  it  may  be  advantageous  to  use  some  of  these  catchment  areas. 

20  In  addition  to  the  above,  a  survey  is  underway  to  collect  information  about  a  network  of 

21  Canadian  watersheds  (CANWANET)  which  have  studies  of  relevance  to  IROWC  (contact 

22  Jacques  Millette).    Most  of  the  studies  are  at  level  2.    See  attached  sample. 

23  6.  ASSOCIATED  RESEARCH  AND  COLLABORATION 

24  Implementation,  testing  and  validation  of  the  proposed  IROWC  methodology  is  best  undertaken 

25  by  building  on,  and  coordinating  with,  related  research  which  is  in  progress  or  planned.    This 

26  section  briefly  identifies  several  associated  research  activities  which  should  contribute  to  further 

27  development  and  application  of  the  IROWC  methodology,  and  also  suggests  a  general  strategy 

28  of  collaboration  among  research  participants. 

29  Research  projects  which  have  potential  association  with  IROWC  at  the  University  of  Guelph 

30  are  listed  below.    This  is  a  partial  and  incomplete  list.    Many  other  research  projects  are  active 

31  at  other  universities  and  research  institutions  across  Canada,  which  should  also  be  considered 

32  for  their  contribution  to  the  methodological  development  of  IROWC.   The  projects  listed  below 
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1  are  selected  and  categorized  relative  to  three  broad  components  of  IROWC  which  require 

2  further  research  and  detailed  data.    The  scale  of  implementation  is  also  shown. 

3  1.     Estimation  of  surplus  nutrients  and  pesticide  residue, 

4  and  their  movement. 

5  Nitrogen  balance  for  agricultural  watersheds  watershed 

6  M.  Goss,  Department  of  Land  Resource  Science 

7  Retention  of  pesticides  by  Ontario  soils  plot 

8  L.J.  Evans,  Department  of  Land  Resource  Science 

9  Movement  of  chemicals  and  bacteria  off/out  of  the  root  zone  field/watershed 

10  D.  Rudolph  Centre  for  Ground  Water  Research  (Waterloo) 

1 1  and  G.  Kachanoski,  Department  of  Land  Resource  Science 

12  (and  other  associated  research) 

13  2.     Contamination  partitioned  by  surface  or  ground  water. 

14  Modelling  the  quality  of  surface  and  tile  drainable  water  field 

15  R.  Rudra,  School  of  Engineering 

16  Field  techniques  for  measuring  the  hydraulic  conductivity  of  soil  field 

17  D.E.  Elrick,  Department  of  Land  Resource  Science 

18  Impact  of  manure  and  fertilizer  on  nitrate  contamination  on  ground  water  field 

19  E.  Beauchamp  and  G,  Kachanoski,  Department  of  Land  Resource  Science 

20  Sustainable  water  use  in  southern  Ontario  regional 

21  R.  Kreutzwiser,  Department  of  Geography 

22  (and  other  associated  research) 

23  3.     Effect  of  management  practices  on  nutrient  balance,  pesticide 

24  fate  and  excess  water. 

25  Effect  of  management  on  unsaturated  transport  properties  of  soil  field 

26  G.  Kachanoski.  Department  of  Land  Resource  Science 

27  Control  of  soil  erosion  and  associated  water  pollution  watershed 

28  T.  Dickenson.  R.  Rudra.  School  of  Engineering 

29  Variable  rate  application  technology  for  N  fertilizers  (effect  on  leaching)  field 

30  G.  Kachanoski.  Department  of  Land  Resource  Science 

31  (and  other  associated  research) 

32  A   strategy   of  collaboration    suggests   that  each   research   participant   contribute    its  area  of 

33  expertise  to  the  development  of  IROWC.    Such  a  strategy  is  briefly  outlined  below. 
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1  1.    Research  and  development    universities,  agricultural  research  stations,  centres  of 

2  of  IROWC:  expertise  (e.g.,  National  Hydrology  Research  Institute) 

3  2.    Collaborative  partnerships:     federal-provincial      agreements      (e.g.,      Green      Plan), 

4  stakeholder       participation       (participatory       research), 

5  international    consultation    (e.g.,    OECD,    Great    Lakes 

6  Commission) 

7  3.     Validation  and  calibration:     water  quality  monitoring  programs 
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2  Figure  1  Nitrogen  Inputs  and  Outputs  for  Wheat  Production  (source:  Juergens-Gschwind 

3  1989) 

4  Figure  2  Climatic  Moisture  Index  -  P-PE  in  the  Prairies  and  Southern  Ontario 

5  Figure  3  IROWC  -  Change  in  Relative  Potential  Contaminant  Concentration  (RPCC)  for 

6  Nitrogen  (1981-  1991)  in  Southern  Ontario 

7  Figure  4a        IROWC  -  Relative  Risk  of  Contamination  in  Southern  Ontario  (Nitrogen  -  1991) 

8  Figure  4b        IROWC  -  Relative  Risk  of  Contamination  in  Southern  Ontario  (Triazine  -  1993) 

9  Figure  4c        IROWC  -  Relative  Risk  of  Contamination  in  Southern  Ontario  (Nitrogen  and 

10  Triazine) 

1 1  Figure  5  Nitrogen  Loading  in  Bamberg  Creek  Watershed  (source:  Ecologistics  Limited 

12  1993T 

13  Figure  6  Nitrogen  Uptake  by  Corn  and  Moisture  Profile,  Elora  Research  Station 

14  Figure  7  Predicted  Annual  Atrazine  Loading  (mg/nr)  Over  Time  at  90  cm  Depth  for  Four 

15  Soils  in  the  Grand  River  Watershed  (source:  de  Jong  et  al.  1994) 

16  Figure  8  IROWC    -    Spatial    Framework    for    Analysis    at    Broad    Level    in    Selected 

17  Watersheds  of  Southern  Ontario 

18  Figure  9  IROWC  -  Relative  Potential  Contaminant  Concentration  (RPCC)  for  Nitrogen 

19  in  Selected  Watersheds  of  Southern  Ontario  1981  -  1991 
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FIG.    6       Nitrogen  Loading  in  Bamberg  Creek  Watershed  (source:  Ecologistics  Limited 
1993) 
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FIG.    7 


Predicted  Annual  Atrazine  Loading  (mg/m2)  Over  Time  at  90  cm  Depth  for  Four 
Soils  in  the  Grand  River  Watershed  (source:  de  Jong  et  al.  1994) 
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CANWANET 
A  network  of  Canadian  watersheds 


Title  of  Study  : 


Wilmot  Watershed  Study 


Name  of  watershed  :    Wilmot  Valley 
Contact  Person  :  Linnell  Edwards 


Location  (province)  :     Prince  Edward  Island 


Objective  of  study  :       Measure  and  model  soil  losses;  develop  soil  conservation  recipes  for  farmers;  develop  a 
local  model;  test  existing  models 


Size  in  hectares  : 
Year  study  started  : 
Description  of  area 

Annual  Cost ; 
Data  collected  ; 


420.00 
04/01/86 


Scale  and  Level 
Present  status  : 


1:1 0000  Level  2 
Ongoing 


Three  sub-watersheds  studied:  140,203  and  420  ha  in  size;  rolling  topography  with 
slopes  up  to  12  %;  annual  rainfall  1040  mm;  soils  are  predominantly  sandy  loams;  crops 
include  potatoes  followed  by  cereals,  tobacco  and  soybeans 


Water  flow,  surface  runoff,  sediment  quality:  N03.  N,  C,  P205,  K20, 


Title  of  Study 


Paired  watershed  study  on  the  Belair  River  Watershed 


Name  of  watershed  :    Belair  River  Watershed 
Contact  Person  :  Jacques  Galichand 


Location  (province)  :     Quebec 


Objective  of  study 


Measure  and  compare  water  quality  of  the  paired  watersheds;  quantify  the  impact  of  the 
treatments  on  water  quality  and  on  the  economy  of  the  untreated  watershed 


Size  in  hectares  :  1 ,030.00 

Year  study  started  :       01/11/94 


Scale  and  Level 
Present  status  : 


1:5000  Level  1  or  2 
Ongoing 


Description  of  area  :     Two  watersheds  studied:  480  (of  which  210  ha  cultivated)  and  550  (of  which  350  ha 

cultivated)  ha  in  size;  topography  mostly  hills,  with  slopes  up  to  10%;  mostly  pastures  and 
hay;  dairy  and  pork  production,  manure  problems:  no  pesticides  used 


Annual  Cost 
Data  collected 


$257,000.00 

Water  flow,  water  temperature,  pH.  air  temperature,  wind  data,  precipitation  including 
snow,  water  quality  parameters  for  samples  taken  weekly  include  N03,  nitrites,  total  N 
NH4-N,  total  P  and  K,  dissolved  P,  BCD,  fecal  coliform  count 


CANWANEI 

A  network  of  Canadian  watersheds 


Title  of  Study : 

Name  of  watershed 
Contact  Person  : 
Objective  of  study  : 

Size  in  hectares  : 
Year  study  started  : 
Description  of  area  : 

Annual  Cost : 
Data  collected  : 


Black  Brook  Experimental  Watershed  Study 


Black  Brook  Experimental  Watershed 


Lien  Chow 


Location  (province)  :     New  Brunswick 


Provide  quantitative  information  on  surface  and  subsurface  runoff;  evaluate  the  impact  of 
various  cropping  systems  environmental  quality;  validate  and  modify  existing  models  for 
N.B. 


1.500.00 
04/01/90 


Scale  and  Level  :  Variable  Levels  1  to  2 

Present  status  :  Ongoing 


Gently  rolling  topography  with  slopes  varying  from  1  to  9%;  mostly  potato  crop,  rotation 
with  grain,  peas  and  hay;  65%  of  area  agricultural  land 


Water  flow  at  8  locations,  surface  runoff,  climatic  data,  water  quality  data:  sediment, 
N-N03,  P,  K,  Mg,  Ca  for  all  locations,  pesticides  at  the  outlet; 


Title  of  Study  : 

Predicting  pesticide  migration  through  soils  of  the  Great  Lakes  Basin 

Name  of  watershed  : 

Grand  River  Watershed 

Contact  Person  : 

Reinder  De  Jong                                        Location  (province)  :    Ontario 

Objective  of  study : 

Develop  methodology  for  predicting;  conduct  preliminary  testing  and  evaluation  of 
methodology 

Size  in  hectares  : 

680,000.00                              Scale  and  Level  :          1 :250000  Level  4 

Year  study  started  : 

04/01/89                                       Present  status :             Ended  in  1994 

Descnption  of  area  : 

Main  crops  are  corn  and  soybeans;  soil  types  range  from  sand  to  silty  clay;  annual 
rainfall  850-975  mm 

Annual  Cost  : 

$100,000  00 

Data  collected  : 

none 
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